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Abstract 

We quantize a generalized electromagnetism in 2 + 1 dimensions which contains a 
higher-order derivative term by using Dirac's method. By introducing auxiliary fields we 
transform the original theory in a lower-order derivative one which can be treated in a 
usual way. 
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I. INTRODUCTION 

In the last years there has been an increasing interest in 2 + 1 dimensions gauge field 
theories [1] in part due to the fact that some of these theories exhibit features which could 
be associated to the high T c superconductors and the fractional quantum Hall effect [2]. 
In such cases the most studied theories have been those which includes a topological mass 
term to the gauge field (Chern-Simons term) that couples to the matter current. 

Besides, the gauge field theories in lower dimensions are a good laboratory to study 
phenomena which are not well understood in four dimensions. 

Particularly, the quantization of Chern-Simons theories has received much attention. 
Dirac's quantization method [3] was used not only in the Abelian massive gauge theory [4] 
but also in the non- Abelian case and in the pure Chern Simons theory [5]. The Abelian 
massive gauge theory was also quantized in a manifestly covariant way [6]. 

On the other hand, theories with higher-order derivatives have also attracted consider- 
able interest. Those terms were already introduced in the past century by Ostrogradskii. In 
the 1940's, aiming to avoid divergences in Maxwell's theory, Podolsky suggested the inclu- 
sion of higher-order derivative terms [7] . More recently, quantization aspects of this theory 
has been analysed by using Dirac's method [8] as well as the Batalin-Fradkin-Vilkovisky 
formalism [9]. This interest rely on the possibility of soften the ultraviolet divergences, 
leading to a possible attenuation of the problem of renormalizability for theories like quan- 
tum gravity. 

Most of the models with higher derivative terms present undesired properties such as 
non-renormalizability and tachyonic massive modes. Despite this, higher-order derivative 
terms have been used in gravitation in order to improve the ultraviolet behavior of the 
Einstein-Hilbert action [10]. Even in supersymmetry and string theory, higher-order terms 
play a certain role. In string theory, for instance, a term proportional to the extrinsic 
curvature of the world sheet was proposed and it has a great influence in the phase structure 
of the theory [11]. The idea was applied to the relativistic particles as well [12]. In 
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supersymmetry, higher-order derivative terms are a useful regularization, that preserves 
supersymmetry. 

In this paper we intend to quantize the Maxwell-Chern-Simons-Podolsky theory by 
using Dirac's method. It was shown that a Chern-Simons term can be generated in Quan- 
tum Electrodynamics in 2 + 1 dimension whenever one integrates the fermionic fields out 
[13-15]. As a matter of fact several other terms can also be generated, including those 
with higher-order derivatives of the gauge field. This can be seen, for example, in the 
derivative expansion method [16], like it was done in reference [14]. In this way, the study 
of properties of the present theory can be helpful to the better understanding of gauge 
effective actions in 2 + 1 dimensions. 

This paper is organized as follows. In section II we define the model and derive the 
propagator. In section III we map the theory into another one which does not present 
higher-order derivative terms and quantize it; in section IV we discuss the results and 
make some remarks. 

II. MAXWELL-CHERN-SIMONS-PODOLSKY THEORY 

We start with the following Lagrangian density 

C = ~ F^FV + °- e^A^d p A u - b l d»F^d x F x » , (1) 

where a, 9 and b are free parameters which permit us taking the appropriate limits. The 
sign of the Podolsky term has been considered in conformity with the original work [7]. 
In order to obtain the propagator one has to add the gauge-fixing Lagrangian 

£ 9 f = ^ (d^f (2) 

to the original one, since (1) is not invertible; where a is the gauge parameter. In such 
case the propagator, in momentum space, is given by 
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G ^ - [( a + 62^2)2^2 _ 9 2 ] { ( a + b2k2 )9^+ 
+a [(a + 6^ 2 ) (a + i + 6 2 ^ 2 - 6 2 ] ^ + ifc^} 



(3) 



which agrees with reference [9] when the parameter 9 goes to zero. It is easy to see that 
the poles of the propagator are defined by the equation 

y 3 + a lV 2 + a 2 y + a 3 = ; (y = k 2 ) , (4) 

where 

2a a? 9 2 

a 1 = — ; a 2 = — ; a 3 = —^ . (5) 

From the above equations we can see that the number of massive poles depends on 
the choice of the parameters. 

By studying the solutions of above equation one sees that, in order to have only real 
roots, it is necessary to impose that the discriminant D = Q 3 + R 2 , where 

« s -I' s 4(t-4 (6) 

be lesser than or equal to zero. So, one have two possibilities: 

i) D = : Here we have three real roots, where two are equal (leading to the appearing 
of ghosts) . So we have that 



mi = 



2R^ 



m = rrin = 



Oi 

3 ' 



Oi 

3 



(7) 
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However, the imposition D = 0, imply that one must have: a 3 = 0, which is a trivial 

3 

solution with 9 = 0; or a 3 = but in this case we get a negative squared mass, revealing 
the existence of tachyonic excitations. 

ii) D < : Now we are faced with three different real roots of Eq.(4). The solutions can 
be written as: 



where 



and 



m. = 



, - l' a \ a l 

2p*cos[-) --, 



cos (f) ± VSsin (|) 



ai 
3 ' 



R 2 + -|a| - a 3 a?/54| 



(8) 



(9a) 



a = arctan 



108y/\a% - a 3 af /54| 



av 



108a 3 



Besides, the imposition D < 0, leaves us with the possibilities: 

< a 3 < §f 



(96) 



1) ai > 0: 

2) ai < 0: g < a 3 < 0. 

In this case however, it is hard to verify if the above values of a 3 and a± correspond 
to positive square masses, at least analitically. Nevertheless, doing a numerical analisys, 
it appears that the following rules hold: 

a) For ai > 0, there is no possible way of keeping the three square masses simultaneously 
positive. 

b) For ai < 0, one observes that along a reasonable range of variation for the parameter 
ai (—0.01 <-> —100), that the region where the three masses are positive, it is always 
greater than that where D < 0. So one can see that the region defined above apparently 
leads to well behaved masses (no tachyons, no ghosts). 
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Unfortunately, the usual case where a = 1, and Podolsky's parameter is negative, 
implies that cii be positive. Consequently the taquyons stays present in the theory. 

III. DIRAC'S QUANTIZATION 

In this section we perform the Dirac quantization through two different ways of gauge 
fixing the theory. 

Looking at the original Lagrangian (1), it is possible to verify that it is equivalently 
described by the Lagrangian density 

C = -- A F^ + e -e^A^d p A v + l -Z^+ b -F^Z^ , (10) 
where Z^ is an auxiliary field and defined by 

Z ltv = d lt Z v -d u Zp , (11) 

is the associated field strength. For sake of curiosity, it is interesting to observe that we 
could also perform a suitable change of variables, capable of disentangle the term coupling 
the field strenghts. Thus we would obtain that 

Z = -\F, U F^ + -l^-(2-sign(-a))Z^ + \z X + 

+ °-e^A^A p - ^—e^Z^Zp - 9 ^e^A^Z p . 
2 4 \a\ va 

In this case it can be seen easily that we would get a Chern-Simons-Maxwell field 
interacting with a Chern-Simons-Proca one through crossed Chern-Simons terms. 

However, we will work with the Lagrangian density (10) instead of that one of equation 
(1). Here we do not have to define the momentum associated with the field time derivative. 
In any case, the introduction of the auxiliary fields duplicates the number of variables we 
are dealing with. 
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The equations of motion are 



a 



+9e^ a F^ = bd^ , 



b d„F^ = Z v 



(12a) 
(126) 



Taking the divergence in equation (12b) and using the antisymmetry property of 
we get 

d v Z v = , (13) 

which is a Lagrangian constraint. 

The equation of motion of the potential has the form 

(a - b 2 n) UA V + 2 6 e^dpAp - (a - 6 2 ["j) d u d»A^ = . (14) 
From the Lagrangian density (10) we have the following primary constraints 



1 = tx q A PS , 

2 = n z w , 

which are, respectively, the momenta associated to the components Aq and Zq. 
The space-components of the momenta are 
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,j 1 - a fi i0 - - 

z 



TV, 



a — — EijAj — b Zio 



7T+ = -b F l0 . 



We can now construct the primary Hamiltonian 



H„ 



d 2 x 

e 



a 

262 



(*i Z ) 2 + 



e 



+ ^ £ ij 7T i Z Aj - EijA diAj 



Z ( ? Z- a „ o 
2 2 4 lJ 



— ZijFij + AiOi + A2^2 



(15a) 
(156) 



(16a) 
(166) 



(17) 



where Aj are the Lagrange multipliers. 

Since the primary constraints must be maintained in the time, their consistency con- 
ditions generates two other constraints (secondary constraints). 



3 = (li = dor* - - EijdiAj « (18a) 



and 

fi 4 = fi 2 = diiTi 2 -Z ^0 . (186) 

The definition of 7r { z in equation (16b) and equation (18b) is nothing more than one 
of the Lagrangian constraints of equation (12b). Furthermore, equation (18b) can be seen 
as a "Gauss law with sources" . 

diEi*-^ , (19) 

where Ei = F i0 . 

It is possible to verify that the consistency of O3 is identically fulfilled, 

3 = . (20) 
On the other hand, consistency of ^4 gives a condition to A2 

4 = -diZ, + A 2 « 

Xi^diZi . (21) 

Up to now we have four constraints. 2 and ^4 are second class constraints and fii and 
Q3 are first class. In order to transform those last into second class constraints we have to do 
a gauge choice. This can be done through two different ways. One can introduce additional 
second class constraints, or break the gauge invariance directly at Lagrangian level by using 
the Faddeev-Popov's trick. Here we will apply the two approaches respectively. 

Equation (14) suggests the following gauge choices 

(a - b 2 U) diAi = , (22) 



diAi = . (23) 
We choose the last one, which is going to be our fifth constraint 

fl 5 = U«0 , (24) 

whose consistency generates a sixth constraint 

Q 6 = O 5 = V 2 A + ^^0 , (25) 

o 

that, in its turn, gives a condition to Ai 

f2 6 = V 2 Ai - ^ « , (26) 

and no more constraints are generated at all. 

It is worth mentioning that the gauge choice (23) can be satisfied if we choose the 
gauge function as [5] 



A=-^d i A i , (27) 



for, starting from the gauge transformation 

A' 3 = A, + djA (28) 

we get 

djA'j = djAj + V 2 A = , (29) 

and, since 

A' = A + d A , (30) 

we have 

V 2 A' = V 2 A -dod t A t = d t F 0t = -^^- , (31) 

o 

where we have used equation (16b). This way we recover the constraint fi 6 , when this is 
strongly imposed, 
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V 2 A>q + = 



(32) 



Now we have only second class constraints and we can invert the constraint matrix. 
The commutation relations among the constraints are 



{n 1 (x),n 6 (y)} = -V 2 S(x-y) , 
{0 2 (x),0 4 (y)} = 8(x - y) , 

{n 3 (x),n 5 ($)} = -v 2 s(x-y) 



(33a) 
(336) 
(33c) 



S(x-y) 



and the others vanish. The inverse of the constraint matrix 

/ v- 2 \_ 

-10 

V" 2 

1 

-V 2 
V-V 2 / 

and the only Dirac brackets among the dynamical variables which do not vanish are 

1 



(34) 



,j - ^ didj)8{x-y) , 



{Z l (x),7i J z (y)} DB = -S i:j 8(x-y) , 
{Zi(x),A (y)} DB = — - ^ diS(x- y) , 
{Zi(f),Z (y)} D _ B _ =d l 5(x-y) 



(35a) 
(356) 
(35c) 
(35d) 



and 



{n l A (x),n J A (y)} D 



B. 



1 



1 



£jfc=2- d k di - e ik —r d k dj 5(x - y) 



From equation (16b) and the Hamilton equation for the momentum tt^ , 



n t z = Z t + bd k F kl , 



(35e) 



(36) 
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we get the "Ampere law in the presence of a source" 

Z 

Ei + e ki d k B = —— , (37) 

where B = e^diAj is the magnetic field. 

Now we perform the Dirac quantization by following the second approach quoted 
before. For this we add the following gauge fixing term in the Lagrangian density (10), 

Cgf = % 2 + a OV), (38) 

where a is an auxiliar field, used to permit that the Lorentz condition be introduced as 
a linear constraint. Its introduction produces the following set of primary second class 
constraints: 

fii = ttq - a ; fi 2 = TTjf ; ^3 = 7r<7) (39) 
whose preservation in time generates the secondary constraints: 



fi 4 = a a - V • A ; fi 5 = V • n z 



fi 6 = a V • An A + %- s ij diAj - V 2 A + • n z , 

£ \0\ 



(40) 



and also eliminate the Lagrange multipliers, 



V 2 A X = - j^V • Z ; A 2 = V • Z ; A 3 = V • ff A + ? e^d z A r (41) 



So we finish with six second class constraints, whose inverse matrix is given by 





( ° 





1 








-V" 2 


\ 
















-2 





1 


-1 








a- 1 










2 








-a- 1 










-2 







2 


















Vv- 2 








o^V" 2 








/ 



5(x-y). (42) 
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After using the constraints strongly, we eliminate the fields A , 7r , ft , ir a , B , ttq. 
Obtaining for the remaining fields the following nonvanishing Dirac brackets: 



B. 



S(x-y); 



1 0*0,- 



• B - 2a|6| V 2 
{Z, , wf } D B = S(x-y); 

{7i l A (x),7i J A (y)} D B = °- (e jk ^2 dkd t -e ik ^ d k d^j S(x-y). 
The reduced Hamiltonian is then written as 



(43a) 
(436) 
(43c) 
(43d) 



H r = ^{X/-A) 2 + a(n A ) 2 + a 6 e^irfAj + ■ n B ) 2 + 
2a 2 

Z 2 CI / _,7,o ft _ . ,_,„•„• CX 6 , o 

- T " 26^* 1 + 2^ + 2H iA) + 

+ |||E iJ Tf^i + \>>\diZjF iJ . (44) 

It is be interesting to note that in this case, as might be expected, the Dirac brackets 
have some dependence in the gauge parameter. Furthermore some of the brackets are 
equal to that obtained using the former approach to fixing the gauge. 

IV. CONCLUSIONS 

We have performed Dirac's quantization of a generalized electromagnetism in 2 + 1 
dimension which contains Maxwell, Chern-Simons and Podolsky's terms. Our results are 
in complete agreement with the literature when the proper limits of the parameters are 
taken. Moreover the freedom to choose these parameters allow us to have distinct massive 
poles for the fields, including the existence of tachyonic modes which, unfortunately, can 
not be eliminated consistently. 

We used the auxiliary field method to reduce the model to a lower-order derivative 
one and consequently pay the price of duplicating the number of dynamical variables. This 
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is an equivalent price one would have to pay if auxiliary field were not used, since one has 
to consider the time-derivative field as well as its conjugated momentum as independent 
variables. 

As a consequence of introducing auxiliary fields we can see the Maxwell equations, 
equations (19) and (37), in a different point of view. The higher-derivative fields can be 
seen as sources to the fields themselves. 

If one considers the effective action to the gauge field generated from the integration 
of fermion fields, the theory treated here should present interesting features. Since Chern- 
Simons and Podolsky's terms can be generated dynamically it would be of interest to 
analyse the dependence with the temperature of Chern-Simons [14,15] and Podolsky's 
parameters, consequently the spectrum of massive excitation. 
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